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Abstract 

We introduce an explicit form of the multi-instanton weight including also instanton- 
anti-instanton interactions for arbitrary in the two-dimensional CP^''~^ model. 
To that end, we use the parametrization of multi-instantons in terms of instanton 
'constituents' which we call 'zindons' for short. We study the statistical mechanics of 
zindons analytically (by means of the Debye-Hiickel approximation) and numerically 
(running a Metropolis algorithm). Though the zindon parametrization allows for a 
complete 'melting' of instantons we find that, through a combination of dynamical 
and purely geometric factors, a dominant portion of topological charge is residing in 
well-separated instantons and anti- instantons. 



1 Introduction 



In the last two decades there have been much evidence that the specific fiuctuations of the 
gluon field carrying topological charge, called instantons, play a very important role in the 
dynamics of QCD, the theory of strong interactions |jl|, 0. Instantons are most probably 
responsible for one of the main features of strong interactions, namely for the spontaneous 
breaking of chiral symmetry ||^. Whether they play a role in the confinement property of 
QCD is, however, not clear. The difficulty in addressing this problem may be connected to 
the fact that, strictly speaking, the multi-instanton solution in the vacuum is not known. 
Instead, in a dilute-gas approximation, one treats the instanton vacuum as a superposi- 
tion of one-instanton solutions. For large distances the one-instanton gluon field falls off 
as Af^{x) ~ which means that correlations between instantons at large distances are 

vanishing in the dilute-gas picture. However, the confinement phenomenon by itself seems 
to indicate that certain correlations of gluon fields do not fall off when the distance is in- 
creased. Therefore, if instantons are of any relevance to confinement, it cannot be seen in 
the dilute-gas approximation: one has to enrich the arsenal of instanton methods. 



This motivates us to study a system where the multi-instanton solution is known, where 
one can find many features of the strong interaction, and which is simple enough to allow for 
analytic solutions in some limits. Such a system is provided by the two-dimensional CP^"^^ 
model. This model is known to possess both asymptotic freedom and instantons. Contrary 
to the d = 4 Yang-Mills theory where the mtt/tz-instanton solutions are not available in an 
explicit form, in the CP^"'^ model they are known explicitly, and for any 'number of colors' 
Nc 0, !§]• Moreover, the multi-instanton weight arising from integrating out quantum oscil- 
lations about instantons is also known @, 0, |^. In addition, the model is exactly solvable 
at large A^'c ^ |^, and the spectrum is known for all Nc |]T^. Therefore, the model is well 



suited to monitor instantons in a controllable way. It should be added that the model has 



been much studied by lattice simulations (see, e.g. |TT|, |T^, |T3[), and instantons have been 
identified there by a 'cooling' procedure [|T^. 

Since the pioneering paper of Fateev, Frolov and Schwartz it has been known that 
the multi-instanton system in the 0(3) sigma model resembles that of a two-dimensional 
Coulomb plasma; it can be exactly bosonized to produce a sine-Gordon theory. In a true vac- 
uum, however, one has both instantons and anti-instantons; such a system can be bosonized 
and solved exactly only at a specific value of the coupling which may not be realistic. 
Therefore, the general case still remains a difficult and unsolved dynamical problem. 



In Sec. we overview the CP^""^ model. In Sec. ^ we derive the general partition function 
for the CP^"^^ in terms of a special form of instanton variables, which are called 'zindons'. In 
Sec. ^ we shall describe an approximation to the partition function called the Debye-Hiickel 
approximation, which is capable to account for many features of the instanton ensemble. 
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However, that approximation is in general insufficient. Therefore, in this paper we combine 
analytical study of the instanton-anti-instanton ensemble with a numerical simulation. In 
Sec. ^we describe the Metropolis algorithm used to study the behavior of the model in Monte 
Carlo simulations. It will turn out that, strictly speaking, the instanton ensemble does not 
exist for Nc = 2 and that a regularization procedure is required to yield meaningful results. 
In Sec. ^ we consider the question whether the system has a thermodynamically stable den- 
sity, i.e., whether the average number of instantons and anti-instantons grows proportional 
to the volume of the system. Finally in Sec. |^ we have a look at the size distribution, both 
seen from the instanton input and from the topological charge density. 



2 The CP^^-^ model 

For a number of colors Nc the main objects of the CP^"'^ model are Nc complex fields 
va{x), a = 1, . . . , Nc- Defining the norm |f p = I]a=i I'^aP one can introduce the normalized 
fields ua, a projector P and the unitary operator G via: 

ua = va/\v\, Pab = uau*b, G = 1-2P, 

p^ = p = p\ G = G\ = 1. (1) 

One can introduce a vector potential and a covariant derivative V^: 

Af, = -{uAdf,u\ - u\df,UA), = 9^ - iA^, {fi = 1, 2) . (2) 

Afj_ is a real field, depending on the x-space variables x and y that can be identified with the 
complex variables z = x + iy and z* = x — iy. The field strength tensor and the topological 
charge density are: 

1 —i 

F^y = df,A^ - d^Af,, q{x) = —e^yF^^ = —e^yd^u*AdvUA- (3) 

47r /TT 

The theory is determined by the action 

S= d'xJ2 \^,ua\' = ^ / ^'a;Tr [d,Gd,G] , (4) 

A=l ^ 

and the partition function is: 

Z = J Dua{x)Du*:^{x)DA^{x)6{\u\^ - 1) exp (^-^ J d^x\V^UA\^^ 

= J DuA{x)Du*A{x)Si\u\^ - 1) exp (-^ J d^xTr[d^Gd,G]\ . (5) 
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This formulates a nonlinear theory of self-interacting fields, where the nonlinearity is forced 
by the condition = 1. There is a topological charge Q in the theory defined as: 

Q = J(fx q{x) = dx^Af, = ^ jd^x e^^F^^ = J d^x e^^ Tr [Gdf,Gd^G] . (6) 
Because of its unitarity G has the property GdG = —dGG, which can be used to show: 

< 1 d^xTi [{d^G ± te^,Gd,G) (d^G ± ie^,Gd,G)] = 16(5 ^ 27iQ) . (7) 

This tells that the minimal action for a given topological charge of a field configuration, i.e. 
a solution of the equation of motion, is obtained if the field satisfies the self-duality equation: 

d^G ± ie^^Gd^G = ^ d^v t itfj^ud^v = 0. (8) 

Introducing the complex derivative = {d^ — idy)/2, it means that one finds two types of 
solutions: 

d^*v = < — > S = 2-kQ (9) 

and 

d^v = < > S = -2nQ , (10) 



which are the Cauchy-Riemann conditions . One calls the first an instanton solution and 
the second one an anti-instanton solution. 



3 'Zindons' 

A general multi-instanton solution of the self-duality Eq. (^) with topological charge Q = 
is a product of monomials in z 0: 

N+ 

vt' = CAY[i^-aM), A = l,...,iV,. (11) 

i=l 

Similarly, a general multi anti-instanton solution of Eq. (|10D with topological charge Q = 
—N- is a product of monomials in the complex conjugate variable z*: 

N- 

vr'=c'^m^*-b\^), A = i,...,N,. (12) 

For a general configuration with instantons and anti-instantons we shall consider a 

N+ Af_ 



product Ansatz |jT5|: 



i=i j=i 
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where aiA and bjA are fixed 2-dim points written as complex numbers. Those points have been 
called 'instanton quarks' or 'instanton constituents' in the past. We suggest a shorter term 
zindon to denote these entities, in analogy with similar objects introduced in the Yang-Mills 
theory fl^. 'Zindon' is Tajik or Persian word for 'castle' or 'prison'; they are of relevance to 
describe confinement in the model. There are, thus, Nc types or 'colors' of the instanton zin- 
dons (denoted by qa) and Nc types of anti-instanton zindons or anti-zindons (denoted by 6^). 

The action and the topological charge for such a multi-instanton-anti-instanton Ansatz are 



S 
2-kQ 



2 / (fx 
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h* ' 
OA] 



AT, 
A=l 



The topological charge can be immediately found from the Cauchy theorem: 



Q = ^ f dx^A^ 

ZTT 



(14) 



(15) 



The geometric meaning of the zindon coordinates becomes especially clear if one takes a 
single-instanton solution, va{z) = ca{z — a a)- Putting it into the action Eq. (|l^) one gets 
the standard form of the instanton profile. 



S = 2ttQ = 2 d'^x 



27r, 



[{x - XqY + p2]2 

where the instanton center xq is given by the center of masses of zindons. 



(16) 



l^A=\ 


ca 


O'Afj, 


l^A= 


1 


CA 


2 



(17) 



while the spread of the field called the instanton size p is given by the spatial dispersion of 
zindons comprising the instanton. 



1^A=\ \^a\ \aA 
E^LiIcaI 



p 





^ _ Ea<b 




2| 




2 


\aA — as 


2 






(Ea 


|ca| 


2^2 



In Eqs.(17, 18) all coordinates xq, a a are understood as 2-dimensional vectors. In the multi- 
instanton case the action and topological charge densities will have extrema not necessarily 
coinciding with the position of the zindon center of masses. Only when a group of zindons of 
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all Nc colors happens to be well spatially isolated from all the rest zindons, in the vicinity of 
this group one can speak about getting a classical instanton profile ([TBI) with the collective 
coordinates given by Eqs.([T^, |T8|), the role of the coefficients ca played by a product of 
separations between a given zindon a^j and all the rest zindons and anti-zindons, 

N+ iV_ 

l^^l^ = n ~ ^Ai'? n - ^Aj?- (19) 

i'^i j=l 

As a matter of fact, it means that there is no need in introducing explicitly extra degrees of 
freedom ca- It should be added that in the thermodynamic limit, with the system volume V 
going to infinity, we expect the number of zindons to be proportional to V. Meanwhile, the 
number of the coefficients ca is fixed (equal to Nc). Therefore, one can neglect this degree 
of freedom as it plays no role in the thermodynamic limit; we shall therefore put all = 1- 

The quantum determinants in the multi-instanton background, defining the instanton weight 
or the measure of integration over collective coordinates qai, has been computed in Ref. 
in case of Nc = 2 (which is the 0(3) sigma model) and in Refs. 0, |^ for arbitrary CP^"'^; 
the results of these references coincide up to notations. 



Da ■ Wa 



N^=2 N+ 

n Ud'aM 

A=l i=l 



nA=i,2n 



i<i' 



n 



case Nc = 


2 




that of the Coulomb 


gas d. 


1 1 2 
\0'Ai — OAj'l 


l(l') 


2 

C{2) 


- a2i' P 



N+ 



(20) 



where A is the renormalization-invariant combination of the bare charge and of the UV 
cutoff M needed to regularize the theory; it is the analog of Aqcd- The numerical coefficient 
C{2) depends on the regularization scheme used; in the Pauli-Villars scheme it has been 
found in Ref. |g to be 

'2\6 3 



C{2) 



- — = 0.04822082 . 



(21) 



We shall absorb this constant in the definition of A. In the one-loop approximation the 
charge in Eq. ( pO]) remains unrenormalized; in two loops it starts to run. It is interesting 
that for accidental reasons (the number of zero modes per unit topological charge, equal 
four, coinciding with the ratio of the second-to-first Gell-Mann-Low coefficients, equal to 
b'/b = ANc/Nc = 4) the coupling g"^ does not 'run' in the two-loop approximation either: 
27i/g'^ in Eq. (pO|) should be just replaced by (6/2)^/^'' = 1. For anti-instanton zindons bAj 
one has a similar expression, with replaced by A^„. We denote the corresponding weight 
by Vb ■ Wb- 
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3.2 Partition function for arbitrary Nc 

The general expression for the muhi-instanton weight Wa in the CP^"^^ model is 0, 



exp 



N, 



^^Inja^i 

i<j A 

/ Nr. 



— - [ (i^x In |f^p(9^ In |t>^| 
Svr J 



Nr 



-^N+ In iY^\cA\']+ -^iV+ + iV+ E In \CA 



\A=1 



A=l 



AT, 



Va 



(22) 



In the following we set \ca\ = 1 for reasons explained above. The second (integral) term in 
Eq. (^) can be analytically computed only for Nc = 2; it gives the Coulomb interaction of 
Eq. (0). For Ai'c > 2 we approximate the integral term by a simpler explicit expression. 
To that end, let us consider a configuration with zindons of Nc colors forming clusters well 
separated from each other, see Fig. |I]. We denote the centers of clusters by x,. 



Xi 



Nr 



(23) 



a=i 



the sizes of the instantons associated with those clusters by pi 



Nc 



Nr 



Xi 



A=1 



^^c A<B 



and the separation between clusters by 

R 



■ij 



(24) 



(25) 



assuming i?^ ^ p? . In such geometry the integral term can be easily evaluated, yielding 



Nc f 

I{aAi} = / (i^a; In |t>^pc}^ In 



Nc 
2 



- N+ N+ 

2j2\nRl + Y.lnp^ + N+ilnNc + l) 

- i<j 



i=l 



+ 



p. 



p 2 



(26) 



We would like to write the functional of zindon positions /{a^j} in such a way that it 

(a) depends only on the separations between individual zindons, a^i — ciBj '■, 

(b) is symmetric under permutations of same-color zindons, {aAi ^ o-Aj); 

(c) reduces to Eq. ( p6|) in the dilute regime; 



7 



(d) at Nc = 2 comes to the exact expression valid for any geometry: 



N+ 

J2 ln(aii-a2,)' + iV+(l-ln2) 



(27) 



The solution to this problem is, naturally, not unique. We shall suggest two forms of ap- 
proximate expressions for J{a^j}, both satisfying the above requirements a-d. The first form 
is 



2 



N 

+ ^iV+(lniV,-l) 



2 2 

Combining it with other terms in ( p2|) we get the following multi-instanton weight: 



(2^ 



exp [3 



N+ r N, 



\n{aAi - aAi'^A^ " X! 

i<i' A=l ^ i,i'=l IA<B 



(29) 



This form of the interactions between zindons belonging to instantons will be used in com- 
puter simulations of the ensemble, see below. This form is, however, inconvenient for analyt- 
ical estimates since it involves many-body interactions represented by lnX]A<_B('^yii — CLBi'Y ■ 
Therefore, we suggest a more simple second form involving only two-body interactions: 



IWai} 



N, 



-—YH H(^A^ - aB,f + ^iV+(lniV, - 1) 



(30) 



The second form reproduces correctly only the leading term in (|2^) in the 'dilute' regime; 
however it is exact for Nf. = 2. Consequently, the multi-instanton weight can be written in 
a compact form: 



Wa = exp 



(3 N, ^+ 

ITTT^T E '^ABHaAi-aBj) A 

^ -L i,j=lA,B=l 



2a2 



(31) 



Here Vab is a projector x A^^^ matrix such that Vab = {Nc — '^)/^c when A = B and 
Tab = ~l/^c when A 7^ 5, so that = V, and V has N^ — l unit eigenvalues and one zero. 



In both forms, (pof ) and (|3lD , we have introduced an 'inverse temperature' /? for future 
convenience, in fact (3=1. We have also borrowed certain powers of the scale A from the 
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overall weight coefficient to make the arguments of the logarithms dimensionless. A com- 
parison of numerical simulations using the weight (|29| ) with analytical calculations using the 
weight (^Tj) shows that the results do not differ significantly; it means that ensembles gener- 
ated with the two weights are to a good accuracy equivalent, hi Nc = 2 they both coincide 
with the exact weight (pOl). For the multi-anti-instanton weight Wh one can use either of 



the forms (poD, with an obvious substitution of the instanton-zindon coordinates aAi by 



anti-zindon coordinates hAi- 

Finally, we turn to the instanton-anti-instanton interaction factor Wah- We assume that 
it is mainly formed by the classical defect of the action, 

Wah = exp (-^Uint{a, b) \ , Uint{a, h) = SAnsatz " 27r(A^+ + A^_), (32) 



where SAnsatz is the action computed on the product Ansatz, (|T^). Again, it is possible to 
evaluate this quantity in a regime, when zindons belonging to definite instantons or anti- 
instantons are grouped together and are well separated from other groups. After some work 
one obtains the classical interaction potential between zindons and anti-zindons in a form of 
'dipole interactions' in two dimensions, valid for any number of colors N^. 

Uint = 871 J2 J[- - Xa)^,ibA - ^h)v^ {^i^v - '^^^ 

all II pairs ^ 

_ ai + ... + aN, _ bi + ... + 6jv, _ _ 

Xa -p^ ; Xfy ^ , iX^ Xljfj,- v'^'^j 



In terms of zindons, this is a many-body interaction, however, Eq. ( p3D can be rewritten 
as a pair interaction of individual zindons, which, naturally, appears to be of a Coulomb 
type. This is one of the great advantages of using the zindon parametrization of instantons. 
We rewrite the dipole interaction (^) as a sum of two-body interactions using the A''^ x 
projector Vab introduced above: 

. N+ iV_ 

f/.„t = -^EE E VABln\aM-bBj\^. (34) 

^^c j=i j=i A,B=1 



c 



In the case when groups of Nc zindons form clusters well separated from groups of 
anti-zindons, this reduces to the dipole interaction (p3|). Notice that the interaction of 
"same-color" zindons is repulsive while that of "different-color" is attractive. At A"c = 2 the 
projector matrix Vab is 

p.-r 4 ( _\ -/ ) . (35) 

so that the interaction becomes just a sum of Coulomb interactions between charges of 
different signs. In this case this form has been known previously |l^ . In Ref . |jl5| it has 
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been shown that this form of the interaction possesses an additional conformal symmetry, 
and one can think that its domain of vahdity in the zindon configuration space is wider than 
the one where it has been actually derived. We shall, thus, use for the mixed instanton- 
anti-instanton weight: 



N+ Af_ AT, 



= exp 2/3A 'Pab In (a^, - 6b,)'A' , A = ^ ■ (36) 



27r 



i=l j=l A,B=1 

The full partition function takes the form (for arbitrary instanton angle 



^-JlJ^jhw^l mi'^A-'— (37) 

with Wab given by (^), and Wa,Wb given either by ( ^9]) or (|3lD . If one uses the form (|3TD 
the interpretation is especially clear: the partition function describes two systems of zindons 
(instanton and anti-instanton ones) experiencing logarithmic Coulomb interactions, whose 
strength is A^^^ — 1 times stronger for same-color zindons than for different- color zindons. 
One has attraction for zindons-anti-zindons of different color and repulsion for zindons- 
anti-zindons of the same color. At iVc = 2 one can think of the ensemble as of that of 
e"*", e~, /i"*", /i^ particles ||T5|]. The interaction of opposite-kind zindons are suppressed by an 
additional factor /5i = 27r/ {g'^Nc). It is actually a running coupling, depending on the scale in 
the problem. Since there is only one scale here, the density of zindons at the thermodynamic 
equilibrium, this /?i should be found self-consistently from the arising equilibrium density at 
the end of the calculations. 



4 The Debye-Hiickel approximation 

To get an insight in the dynamics of the system by simple analytical methods we consider 
here the so-called Debye-Hiickel approximation to the partition function (^). Let us consider 
first a general case of a statistical ensemble of F kinds of particles with two-body interactions 
given by: 

Uint = J2 '^mtiXi - Xj) . (38) 

The number of particles of kind / is Nf, f = 1, . . . , F. One writes the partition function for 
fixed numbers of particles as: 

^ = = / n n . (39) 
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where T is the Helmholtz free energy and are coordinates of particles of kind /. The 
one-particle density of particles of kind / is defined as: 



11/- 



dx(nf{x() = 1 , (40) 



f dxi... dx{^ dx{dxl ■ ■ ■ dx%^ n n dx'l^e-'^^"^' . (41) 



h?^/ j. =1 



which differs from the full partition function by that one does not integrate over the coor- 
dinates of one particle of kind /. Similarly, one can introduce the two-particle density by 
avoiding integrating over coordinates of particle / and of particle g: 

M4,xf) = -^^,___^^, 

This function can be written as 

rifM^^) = nf{x()ng{x',)e-^^fM-^'^ , (42) 

which serves as a definition of the two-particle correlation function ujfg{xi,X2)- If this func- 
tion is zero the particles are not correlated. Introducing, in the similar fashion, higher-order 
densities, one can derive an (infinite) chain of equation for those densities. The so-called 
correlation energy is: 



Ecorr = E / / dx^dx2Ufgixi - X2)e-^'' ^ ^^"'^ . (43) 

fg 



Knowing the correlation energy as function of temperature one can restore the Helmholtz 
free energy from the general relation: 

i^corr = . (44) 

The Debye-Hiickel approximation consists in assuming that Pwfg is small and linearizing 
the equations [19|. It is usually justified, when the temperature is high, /5 — > 0. In the 
Debye-Hiickel approximation it is also convenient to pass to the Fourier transforms of the 
pair interaction potential ujg and of the correlation function Ufg. The basic equations from 
where one finds the correlation functions take the form: 



^fg(k) = Ufgik) + Pj2^l^ghi-k)uJfhik) + Ufh{k)ujghi-k)] . (45) 



In our case the particles are zindons belonging either to instantons or to anti-instantons. 
Let us denote uaiA2,^AiA2 the interaction and the correlation of zindons of color Ai, 
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belonging to instantons. Similarly, Ub-^b21^BiB2 refer to zindons belonging to anti-instantons. 
and u^Bi^AB refer to instanton zindons of color A and anti-instanton zindons of color B. 
The approximation is applicable only if all interactions are two-body. For that reason we 
choose to work with the instanton weights given by Eqs. (|3TD , (|36|) : 



^ - ^E_^(^/f«f''A-'"-"-'~. (46) 

Since our interactions are all of the logarithmic type the Fourier transforms of u^g can be 
written as: 

47r 1 

W^iAal^) = ^AiAap-, ^^^1^1 = !, ^AiAa = ~ ^ _^ if Ai ^ A2, 

47r 1 - - 

Ub.bM = VB,B2-j^, ^^BiBi = l, Vb^B2 = - p^ _i if 517^^2, 

Air N —1 1 

UAsik) = Vabj^, ^aa = 2A^^, t;^^ = -2/3i— if A ^ B . (47) 

It is now convenient to change the notation for the u correlation functions to: 

UJaA = <-^l, ^ AA' = Cl'2, '-^BB = ^3, ^ B& = UJ4, UJj^B = ^5, ^ A& = UJq- (48) 
AytA' B^B' B'^A 

The above relations for the w's induce similar relations for the cu's: 

^2 = -^^1, = -J^^-^s, = -]^^5 . (49) 

With the new notations we get the following set of linear equations for the three independent 
correlation functions cui 2,3: 

/3uJi = -7[1 + n+(3uji + /3'h_l3u5\ , 

Puj5 = --f[p' + ^P'P{h+Ui + n^LJ3) + ^Puj5{n+ + n.)], (50) 



where 



n I = = n+ , = = . 51 

We get the following solutions for the set of linear equations: 

-47r[A;2 + 47r/3n_(l -/3'2)] 



(A;2 + fi:f)(A;2 + /t2) 
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(F + Kf)(A;2 + /s:2) 



2ttI3 



(52) 



Returning to the x-space we observe that all correlation functions are exponentially de- 
creasing at large separations between zindons, as contrasted to the original logarithmically 
growing interactions. This phenomenon is usually referred to as 'Debye screening'; in our 
case it takes the particular form of Eq. (|52D. 



Using Eqs. ( ^ , ^1D one can reconstruct the partition function for a fixed number of instantons 
and anti-instantons: 



Z{N+,N. 



exp 



1 



4(iV, 



-V 



In 



A2 



477/3. /n+n_ (1 



+ 1 



— -\/(n+ — n_Y + Ah+h_3''^ In — 



{VA 



2\NciN++N_) 



A = 2Ae-^^ . 

At N_^=N_ = N/2 corresponding to the instanton angle 6* = it simplifies to: 



(53) 



Z{N) 



exp 



1 - In 



MNc-i)v, 



{VA 



2\NcN 



m 



2Nc 



(54) 



Finally, one can find the optimal density of instantons and anti-instantons by maximizing 
the above partition function in A^. For a fixed (3' and A we get then the following density of 
instantons plus anti-instantons: 



n = {n_ + n+) 



<N> 



(3')— An (3 



A^. 



2 ^^l 



fJN^ 



{iVc~l)2-/3JV^ 



(55) 

Notice that the density is stable in the thermodynamic limit V ^ oo. At \j3'\ = 1 corre- 
sponding to the value of the instanton-anti- instanton interaction constant j3i = =l= 2(jv^-i) 
system becomes unstable. This could be anticipated from direct inspection of (^): at that 
value of jSi the attraction between zindons of different color belonging to instantons and to 
anti-instantons exceeds the repulsion between same-kind zindons, and the system collapses. 
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It should be stressed that all quantities involved (like the instanton density n, the screening 
masses ki^2 etc.) are stable in the limit Nc — >■ oo. This should be contrasted to the old con- 
jecture by Witten [§] that instantons die out in the large A^^^ limit. Though this conjecture 
has been criticized later on 



20, 21 1 it remains a widely-believed prejudice. 



5 The Metropolis algorithm and regularization proce- 
dure 

The numerical treatment of instantons in the CP^"'^ model is problematic for N^. = 2, where 
the theory is equal to the two-dimensional 0(3) model. The partition function in this case 
involves a product of integrals of the type / dX^/X^ and is therefore ultraviolet divergent. 
There are two ways of imposing a regularization. The first consists in decreasing the powers 
in the interactions (^). Physically, this is equivalent to the introduction of a temperature 
T = 1/(3, see Eqs. ([29|) , (pT|) , and (^). We reach the physical limit by taking (3 —>■ 1. For 
small values of /5 we expect the model to be consistent with the predictions of the Debye- 
Hiickel approximation. The second possibility is using an ultraviolet cutoff e. Technically, 
this can be easiest done by changing the logarithmic interaction between zindons: 

ln{aAi - aAi'Y In (a^i - aAi'Y + ^ , etc . (56) 

In both cases the physical limit T ^ 1 and e does not exist for the free energy, however 
one can get meaningful results for physical quantities that are connected to derivatives of the 
free energy. An unpleasant point is that this e regularization destroys the scaling properties 
of the theory. We can restore the scaling property by using e' = e/L clS cl gularization 
parameter. Then we restore the original scaling property which we express here in terms of 
the free energy F: 

F(e', AT, L) = In Z{e', N,L) = N In(L^) + In z{e, N) . (57) 

Note that the reduced partition function z{e', N) does no longer depend on L. It is exactly 
this scaling property which will be important later when we try to reconstruct the average 
density of instantons and anti-instantons in the system. We should emphasize that the price 
we have to pay for the restoration of the scaling behavior is that the two limits e' — > 
and L — »• oo cannot be interchanged, as with growing L the effective cutoff e = e'L becomes 
larger than the average separation between zindons, (R) = 1/ ^/n. However, when e' becomes 
smaller the physical region where e = e'L < (R) covers larger portions of the configuration 
space. 

To study the behavior of the system numerically we first switch off the interaction between 
instantons and anti-instantons i.e. we put /3i = 0. We then have for Nc = 2 the case of 
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a simple Coulomb gas, where the two different colors interact in the same way as the two 
different electric charges. In the case of larger Nc we have a more complicated 'multicolor' 
interaction which would then require a picture with Nc different kind of charges. 

The variables of integrations for the theory are the zindon coordinates aAi and A 
direct numerical calculation of the partition function (^) is technically not possible, be- 
cause the integrand has very sharp maxima: a small variation of the zindon coordinates 
changes its value by orders of magnitude. Instead of this one can apply a Metropolis algo- 
rithm: Starting with a random configuration of zindons, step by step the position of one 
zindon is varied by a small amount. If the new configuration has a larger value of the 
weight WaWbWab then this configuration is accepted, if not then it is only accepted with the 
probability WaWbWab{^G'Vi')/waWbWab{o\d). As Fig. ^ shows, the plateau of important config- 
urations is reached only after 250.000 Metropolis sweeps for N^. = 3,4. In the case = 2 
the number of steps necessary to reach the plateau depends on the regularization parameter 
e': It is large when e' is small. In Fig. |^ we have chosen e'^ = 0.001. Already after 20.000 
steps one is very close to the plateau here. It is now interesting to make a snapshot of the 
distribution of the zindons at the plateau region. In Fig. ^ we show snapshots of the zindon 
ensemble at the plateau region after 500.000 Metropolis steps with j3i = 0. Each symbol 
represents zindons of one color. One observes clearly, that for = 2 zindons of different 
colors tend to condense into neutral pairs that can be very densely packed. This is the 
reason for the fact that the partition function is divergent and needs to be regularized in 
the way described above. Such a condensation does not take place for higher colors Nc = 3,4. 

We can put this in a more quantitative way by plotting the free energy F = InZ as a 
function of the regularization parameter e, see Fig. ^. As it should be expected by simple 
power counting the free energy for the case Nc = 2 diverges proportional to ln(ln(l/e)), 
while it is finite in the case of higher Nc. In principle such a divergence does not pose a real 
problem since physical quantities are only related to logarithmic derivatives of the partition 
function, but it does affect the effectiveness of the numerical studies considerably. In fact, 
with the computer powers used, it was not possible to go to e'^ smaller than 0.001, which 
is not very small. The smaller e'^ becomes the larger are the fluctuations in the partition 
functions until a point is reached where we do not get any equilibrium within a reasonable 
time at all. We have to conclude from this that the Nc = 2 case is not a good system to be 
studied numerically. 



6 The stability of the instanton— anti-instanton ensem- 
ble 
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6.1 The case iV^ = 2 



The first question is whether this system has a thermodynamic equilibrium, i.e. whether its 
density does not change with the system volume L^, when L is large. We start again with the 
Coulomb gas of two noninteracting ensembles of instantons and anti-instantons, i.e. {(3i = 0) 
and Nc = 2. Here it is sufficient to consider the instantons alone. In this case we simply 
have = and we plot the logarithm of the partition function, i.e. the free energy as a 
function of the number of instantons A^. The ideal situation would be now that the number 
-^max; where the partition function has a maximum, is proportional to L^. Then one could 
continue this situation to an infinite instanton ensemble with the same density. In Fig. § we 
show the partition function for the Coulomb gas as a function of the number of instantons 
N for various box lengths L. The simulation is done with the UV cutoff e'^ = 0.001. One 
sees that the free energy is plagued by large fluctuations that correspond to fluctuations in 
orders of magnitude in the partition function. A direct determination of the maximum of 
the free energy is not reliable. Instead, one can exploit the expected scaling behavior. If the 
system has a stable density then the free energy must be of the form: 

F(e', N, L) = In Z(e', A^, L) = N{\n{L'^) - [\n{N) - 1 - ln(r2+(e'))]} + a , (58) 

with the two fit constants a and n+{e'), and e' = e/L being fixed in order to preserve the 
scaling behavior of the initially unregularized theory. Eq. ( pSf ) is obtained in the following 
way: From the inspection of Eq. (^) it is clear that the free energy depends on L via 
A^ln(L^). Requiring now that the number A^ at the maximum of the free energy F is 
proportional to leads to a differential equation with the solution given by Eq. (^8l). 
Eq. (|58|) has the advantage that n^{e') gives immediately the density of the system. The 
symbol n+ refers to the instanton density. Unfortunately, the system does not exhibit a 
stable density. The reason is, again, the tendency of the system to condense into neutral 
pairs corresponding to small-size instantons. When it happens, it is not reasonable anymore 
to write down the particle-identity factor 1/(A^!)^ in the partition function: it should be 
rather replaced by the first power of the factorial. Therefore, we allow for an additional 
'conbinatorical defect' term in the fit, which reads alnA^!, so that the full fit function has 
now the form: 

Ffit(e',Ar,L) = N{\n{L'^) - [\n{N) - 1 - ln(r2+(e'))]} + aln A^! + a . (59) 

The smooth solid lines in Fig. |^ are the fitted curves, which fit quite well to the heavily 
fluctuating free energy. One should expect that the defect becomes smaller when e' becomes 
smaller. As shown in Tab. |I|, one sees indeed such a behavior. 

Now arises the question how reliable our simulations are, because so far we did not find 
strictly speaking a stable system. The situation is also not very comfortable because in the 
limit e' = 0, where a stable phase should exist the density n+ is divergent for Nf. = 2. We 
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therefore study now a different regularization method: It consists in introducing a tempera- 
ture T = in the zindon interactions, with /3 < 1. We reach the physical hmit by taking 
/3 — >■ 1. For small values of (3 we expect behavior of the system to be consistent with the 
predictions of the Debye-Hiickel approximation. In Fig. ^ we show the free energy for /3 = 0.1 
versus the number of instantons A^. The solid lines are the fits using the form: 

-Pfit, Debye 

(iV, L) = iV{ln(L2) - [In(iV) - 1 - ln(n+,Debye)]} + a . (60) 

We find a stable density for various box lengths L. The five curves yield a value for 
n-i- = 1.06787 ± 0.0024, which is in agreement with the prediction of the Debye-Hiickel 
approximation: (p5D gives Debye = 1-0370. 



We now switch on the interaction between instantons and anti-instantons putting /5i > 0. 
Again we can only determine a kind of stable density n+(e'), if we allow for a defect param- 
eter. Tab. 1^ shows the results for n+(e') for three different values of e' as a function of Pi. 
In case of /?i = we can compare the density to the one in Tab. [l|, both are results of inde- 
pendent simulations. A comparison of the upper line in Tab. ^ with corresponding entries 
in Tab. ^ gives an idea of the statistical error; it is not enormous, but certainly not negligible. 

We observe that the density rises only slowly with This is in agreement with the results 
of the Debye-Hiickel approximation, which also predicts a slow rising of the density with jSi. 
For Nc = 2, the equilibrium density is given by (see Eq. (0)): 



/(/5i)7r/5e2^-2f 
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l+|3l , ^ , l-fa 



/(A) = (i + A)— (i-/3i)— . (61) 



The function /(/5i) varies between 1 and 2 for /3i varying between and 1. For the real CP^ 
model, when /3 = 1, we can expect at least a behavior linear in ln/(/3i) for the logarithm 
of the density lnn+. Fig. shows, indeed, that we can obtain a good linear fit, using this 
assumption. The numerical results for the linear fit can be read off from Tab. |^. The growth 
of the intercept a when e'^ goes to zero reflects the fact that in this theory the density is 
divergent. The growth of the slope h is also natural, since at e' — *• we expect h = 1 from 
Eq. (|6T| ) . In general one sees that the effect of the instanton-anti-instanton coupling on the 
density is comparatively small as predicted also by the Debye-Hiickel approximation. 



Finally, let us discuss the running of the instanton-anti-instanton coupling From its 
derivation in Sec. it is seen that j3i = vr/gf^ ( at Nc = 2) is related to the bare coupling 
constant of the theory g"^. However, physical quantities cannot depend on the bare couplings, 
but only on the renormalized ones. It means that by considering quantum corrections to the 
classical instanton-anti-instanton interactions one would be able to observe that the bare 
coupling constant g'^ (or j3i) is replaced by the running coupling depending on the charac- 
teristic scale in the problem at hand. In our case the only dimensional scale is the average 
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separation between instantons and anti-instantons or, else, the density of zindons, n. From 
renormalization-group arguments one can write the 1-loop formula 

/3i = ^(ln^ + c), (62) 

however, the constant C remains unknown. In order to determine C one has to perform the 
renormalization of the instanton-anti-instantons interaction with a very high precision; this 
has not been done. 

Keeping in mind that the dependence of the instanton density on f3i is very mild in the 
whole range < /?i < 1 (see Fig. ^ and Tab. 0) and that the dependence of Pi on the density 
is also quite weak (see (|62D), and given the theoretical uncertainty in the constant C, we 
cannot reliably determine the density self-consistently by taking into account the running 
of In the case of Nc = 2 by far a more important factor affecting the density is the 
necessary ultraviolet cutoff e. 

6.2 The case iV^ = 3, 4 

In the case Nc = 3,4 we are not plagued with the necessity to introduce a regulator. However, 
we are confronted now with enormous fluctuations. In Fig. ^ the free energy is plotted for 
several values of f3i. Already in the free energy the fluctuations are of an order of magnitude, 
so that a reasonable determination of the equilibrium density is not possible. However, we 
could at least try to compare the situation to the one in the Debye-Hiickel approximation 
for Nc = 3, 4. First of all it is necessary to mention that our partition function ( ^Tj) is not 
based on two-body interactions only, and is therefore not directly comparable to the one used 
to derive the density (pSj). But, as we will see now, some of the features are nevertheless 
strikingly in parallel. First of all it turns out that for all values of f3i not too close to the 
critical value, the free energy does not change much (as compared to the extremely large 
fluctuations). This is in parallel to the fact that the equilibrium density given by the Debye- 
Hiickel approximation grows only slowly with f3i, see Fig. ^. Another feature of the Debye- 
Hiickel result for the equilibrium density is that it is comparatively large for Nc = 2 and then 
falls rapidly to the asymptotic value n^{Nc = oo, /? = 1, /?i = 0) = n^^'^ exp(27£;/7) ~ 1.3888. 
The second point, where there is something in parallel between the Debye-Hiickel result and 
the simulations, concerns the critical behavior. The Debye-Hiickel formula for the density 
becomes meaningless when jSi > Nc/2{Nc — 1). Hence for A^^c = 3 we find jScru = 3/4 and for 
A'c = 4 we find /3crit = 2/3. In the simulations it is seen (Fig. |^) that for values of /3i ~ 1 
the system behaves abnormally and shows clear signs of a collapse. 
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7 Size distribution 



A very important question is what is the size distribution of instantons in the ensemble. 
Is the average size (p) smaller or larger than the average separation between instantons 
(R) = n~^/^? In the former case one can say that instantons are dilute while in the latter 
case one says that instantons 'melt', so that individual instantons have not much sense. 

We are now in a position to study this problem quantitatively. The zindon parametriza- 
tion of instantons is an ideal tool for that since it allows for both extreme cases. If zindons 
of different 'color' are uniformly distributed in space that would mean the instantons 'melt'. 
If zindons of Nc colors tend to form well-isolated color-neutral clusters, it means instantons 
are dilute. The size distribution of instantons has been measured on the lattice for = 2 
TM after a cooling procedure has been applied to remove short-range fluctuations. (For a 



recent comparison of different cooling techniques on the CP^''^^ model see also |2^.) A 
very delicate question in connection to this is, how to identify instantons on the lattice. 
Several procedures are in use In this paper we want to compare two ways of extracting 
the instanton content. The first one, which we call the 'geometric' method, makes use of 
the zindon picture. From the sample of all zindons the color neutral iVc-plet which has the 
least dispersion about the common center is chosen to be an instanton. Then all zindons 
belonging to this instanton are removed from the sample, and the procedure is iterated to 
find the next group of zindons that has the smallest dispersion and so on, until the whole 
sample is grouped into instantons. The second procedure, which we call the 'lattice' method 
as it is more close to the ones used in lattice studies, is looking for local maxima of the 
topological charge density g(x) on a N^^^^ x iVgrid lattice. An instanton is assumed if a grid 
point has a larger topological charge than its surrounding 8 nearest neighbors. In order to 
improve the accuracy in finding the center points and the sizes, we construct an interpolating 
polynomial, 

g(x, y) = ax^y^ + hx^y + cxy"^ + dx'^ + ey"^ + fxy + gx + hy + j, (63) 

calculated from the 9 grid points consisting of the grid point under consideration and its 8 
neighbors. From this interpolating function we find coordinates of the local extrema xe- An 
extremum is accepted if it lies inside the 3x3 grid used for the interpolation. In order to 
classify whether we have found a local maximum, minimum or a saddle point we construct 
the matrix of second derivatives at the extremum xe'- 



\ dx J dx dy 

d_d_ (d_\^ 
\ dx dy \dy ) 



q{xE) . (64) 



If the two eigenvalues Ai and A2 of this matrix have a different sign then is a saddle 
point and is rejected. For negative Ai and A2 we have a local maximum and respectively a 
local minimum if both Ai and A2 are positive. The two eigenvalues give at the same time a 



19 



measure for the size pf = \J^/\2'n\i\, i = 1,2 (they should coincide for an ideal spherically- 
symmetric instanton). This size can be compared to the one derived from the value of the 
topological charge at the maximum = 2/|27rg(x£;)|. Only if all three sizes are smaller 
than L/2 is the instantons or anti- instanton accepted and the size is then taken to be the 
geometrical mean of the three definitions, p = {pip2PzY^^ ■ It should be noted that usually 
Pi, 2,3 are not drastically different. We should also stress that a reliable identification of all 
local extrema is achieved only with a fine grid. With a 100x100 grid we find almost the 
same number of maxima (instantons) and minima (anti-instantons) as we put in 'by hands' 
by fixing the number of zindons and anti-zindons. 



We first test both methods, the 'geometric' and the 'lattice', in a toy model where the 
positions of zindons and anti-zindons are random: the interactions are switched off com- 
pletely (/3 = and (3i = 0) so that the distribution of zindons in space is given by a 
fiat measure JJcPaAiYld'^bAj- We generate configurations of zindons with this measure and 
compute the topological charge density q{x) of a given configuration from the integrand of 
Eq. (|T¥). Instantons and anti-instantons are then identified by the two methods described 
above. Taking many configurations we compute the distribution of the instanton sizes, see 
Fig. M. 



For the fiat measure the small-p side of the distribution can be immediatelly evaluated 
from dimensions to be i/(p) ~ p^^"^^. The geometric distribution (shown by histograms) 
follows this prediction whereas the lattice method strongly deviates from it. [One can mimic 
the smooth geometric distribution by the lattice method by taking a crude grid but then 
many instantons are lost as they "fall through" the grid.] The fact is that the topological 
density computed from randomly distributed zindons shows many sharp and narrow peaks. 
The larger number of instantons we take (at fixed density) the more sharp the peaks become: 
it clearly indicates a non-thermodynamic behavior of the system. 



This pathological behavior is due to the utmost instability of the product Ansatz ([13|) in the 
case when the interactions are swithched off. Indeed, the sizes of instantons are determined 
by Eq. (0) where the role of the weights ca is played by long products (pJ]). These products 
are extremely unstable: sufficient to change the positions of very distant zindons, and cas 
will change. With zindons distributed according to the fiat measure c^'s fluctuate by orders 
of magnitude, and the more zindons are taken the stronger they fluctuate. It means that 
there will always be a zindon of a certain color whose weight ca is much much larger than 
those of the other zindons forming an instanton. According to Eqs. (|17| 0) the extremum of 
the topological density will then coincide with the position of the zindon of that particular 
color, and the size of the corresponding instanton will be close to zero. This is why we see 
too many small-size instantons by the lattice method. See also the discussion around Fig. ^ 
below. 
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When we switch on the interactions zindons get correlated and the product Ansatz should 
not be pathological any more. In particular, from the experience of the Debye-Hiickel ap- 
proximation we know that the correlation functions of zindons get exponentially suppressed 
at large separations. Therefore, we expect that the topological density will not have un- 
physical sharp peaks at small p as in the case of the flat measure, and the two methods of 
instanton identification should come closer. In practice, however, one cannot be sure before- 
hand that the thermodynamic equilibrium is fully reached, given the limited computer time 
for configuration generation. 

We switch on the interactions, first, only of zindons belonging to same- kind instantons 
{(3i = 0) and then adding interactions between instantons and anti-instantons = 0.5), 
see Fig. [Tl|. We show only the cases A^c = 3,4 as for Nc = 2 the size distribution is cutoff 
dependent. For the 'lattice' method we use again the 100 x 100 grid. Unfortunately, the 
fluctuations of the partition function in the simulations are so large that we have to apply a 
trick in order to get meaningful results. We let the Metropolis algorithm work until a plateau 
is reached. Then we make many measurements of the size distribution at the plateau, all 
weighted with a weight averaged over the plateau. We repeat this for several starts but we 
essentially see the size distribution for the plateau possessing the largest free energy. 

The price we have to pay for this somewhat truncated simulation is that the unphysical 
peaks at small p, as determined from the lattice method, are still present though they are 
not so awkward as in the non-interacting case. As Nc increases the small-p peak shrinks, 
therefore one can speculate that at large A^ it is easier to reach the thermodynamic equilib- 
rium. Except for very small p the two methods of instanton identification seem to produce 
rather similar distributions. 

The geometric size distribution exhibits the expected behavior at all p. In particular, at 
small p it behaves as z/(p) ~ p^~^ where b = Nc is the 1-loop Gell-Mann-Low coefficient, 
which is to be expected from dimensional analysis. 

We observe that the instanton-anti-instanton interaction does not infiuence the size dis- 
tributions significantly. This is seen both from geometric and lattice methods by comparing 
curves in Fig. [ll] with (3i = and Pi = 0.5. The case Pi = 0.5 has a slightly broader 
distribution. Qualitatively, it should be anticipated since the interactions between zindons 
and anti-zindons tend to 'pull out' individual zindons out of their neutral clusters, however, 
the effect is small. 

As to the interactions of the same-kind zindons, it clearly leads to smaller instanton sizes; this 
could be also anticipated. Indeed, the zindon interactions are such that same-color zindons 
are repulsive while different-color zindons are, on the average, attractive. Such interactions 
lead to preferred clustering of zindons into color-neutral groups. (This is especially clear 
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for Nc = 2 where the Coulomb interaction leads, as we have seen, to a collapse of zindons 
into neutral pairs as one removes the ultraviolet cutoff, e — 0). In Fig. 11 we show by 
arrows at the top of the plots the positions of the 'geometric' distributions maxima in the 
non-interacting case of Fig. |l^: both for A^^^ = 3 and 4 the maxima are noticeably shifted to 
smaller sizes when one switches on the interactions. The 'lattice' method shows this shift 
only in the case iVc = 4 where the unphysical small-size peaks of the topological density are 
suppressed. 



In Fig. |T2| we plot the topological density obtained from a typical configuration of zindons 
at Nc = 3 with all interactions taken into account. We see that the peaks of the topological 
density in general correlate with the positions of the centers of instantons found by the 'ge- 
ometric' method (shown by triangles in Fig. 0). In some cases, however, the peaks of the 
topological charge density sit not in the centers of the 3-plet of zindons forming an instanton 
but rather near one of the three zindons. This is the consequence of the instability of the 
weight coefficients ca, discussed above, and is the actual cause of the small-p peaks in the 
'lattice' size distributions. 



Finally, we would like to draw attention to the fact that in Fig. |n| we have plotted the instan- 
ton sizes in units of the average separation between instantons, defined as (R) = {N/V)~^^'^ . 
It is remarkable that both for A^^^ = 3 and 4 and independently of the method used to identify 
instantons, the average instanton size is definitely smaller than the average separation. It 
means that instantons formed of A'c-plets of zindons are, on the average, dilute. Though one 
can always find a huge instanton overlapping with many others, the probability of such an 
event is relatively low. 



8 Conclusions 

We have formulated the statistical mechanics of instantons and anti-instantons in the d = 2 
CP^"'^ models in terms of their 'constituents' which we call 'zindons'. We have derived the 
interactions of same-kind and opposite-kind zindons for arbitrary Nc- At Nc = 2 they come 
to logarithmic Coulomb interactions known previously. At A^c > 2 the interactions are of 
a many-body type though one can reasonably approximate them by two body logarithmic 
interactions with charges depending on the 'color' of zindons. 

At high temperatures the system can be studied analytically by ways of the Debye-Hiickel 
approximation, at any Nc. The physical case corresponds to the temperature P = I and can 
be studied by a Metropolis-type simulation. The numerical study shows that the Debye- 
Hiickel approximation is, as expected, accurate at small /?, but is qualitatively valid even 
at (3 = 1. Both analytical an numerical studies indicate that the effect of instanton-anti- 
instanton interactions (in our language the zindon-anti-zindon interactions) is not significant: 
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the interaction of same-kind zindons or, else, the nmlti-instanton measure, is by far more 
important for the statistical mechanics of the instantons-anti-instanton ensemble. 

At A^c = 2 the system is, strictly speaking, collapsing into instantons of zero size: to reg- 
ularize its behavior one needs to introduce an ultraviolet cutoff. We have studied in detail 
the statistical mechanics with varying cutoff parameter: the numerical simulations are fully 
compatible with the expectations. 

At Ac = 3, 4 the system does thermodynamically exist with no cutoff, however, it exhibits 
enormous fluctuations which undermine accurate measurements of physical quantities. Nev- 
ertheless, we have managed to study the instanton size distribution introducing two methods 
of its determination, 'geometric' and 'lattice'. Using the 'lattice' method with a very fine 
grid we observe many unphysical small size instantons for Ac = 3, which stem from the 
instabihty of the product Ansatz. Though this effect is believed to disappear when the sys- 
tem approaches the true statistical equilibrium it is very difficult to get rid of it in practice 
because of gigantic fluctuations in the system. The effect is even more pronounced for the 
flat measure but tends to vanish when Ac is increased. 

Though the zindon parametrization of instantons and of their interactions allows for com- 
plete 'melting' of instantons and is quite opposite in spirit to dilute gas Ansatze, we observe 
that zindons, nevertheless, tend to form 'color-neutral' clusters which can be identified with 
well-isolated instantons. This effect is due to a combination of two different factors both 
supporting clustering. One factor is the interactions: same-color zindons are strongly repul- 
sive while different- color zindons are attractive. The second factor is pure geometry: even 
with a flat measure the probability to combine A^c zindons into a neutral cluster smaller than 
the average separation is quite sizeable. Both these factors are expected to be even stronger 
in four dimensions appropriate for the Yang-Mills instantons. 
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density n+ 


average distance 




N 

'max 


defect a 


0.1 


5.095 ± 0.021 


0.443 


1.401 


10 


0.839 


0.03 


12.707 ± 0.107 


0.281 


1.620 


33 


0.777 


0.01 


23.501 ± 0.211 


0.206 


2.063 


100 


0.672 


0.003 


42.723 ± 0.673 


0.153 


2.793 


333 


0.571 


0.001 


64.117 ± 1.193 


0.125 


3.949 


1000 


0.481 


0.0003 


135.207 ± 4.299 


0.086 


4.965 


3333 


0.518 



Table 1: Dependence of the zindon density in the Coulomb gas on the cutoff parameter. 
The error is given by the standard deviation of the n+'s obtained by separate fits to five 
measured curves. The value of n+ itself is the value one obtains if one performs a combined 
fit with one n+ to all the measured curves. The real value of n+, however, is determined by 
the error in the defect a, to gain experience of its influence one may compare to Tab. |. 
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I3i 


e' = 0.1 


e'^ = 0.01 


e'^ = 0.001 




a = 0.7539 


a = .5785 


a = 0.45992 




[A = -0.036) 


[A = -0.035) 


A = -0.013 


0.0 


6.946 


26.314 


74.948 


0.1 


7.062 


27.614 


75.371 


0.2 


6.945 


27.650 


75.250 


0.3 


7.060 


28.229 


74.075 


0.4 


7.080 


27.711 


81.660 


0.5 


7.164 


27.707 


77.512 


0.6 


7.006 


29.171 


80.946 


0.7 


7.149 


28.706 


82.323 


0.8 


7.306 


30.013 


91.595 


0.9 


7.395 


31.705 


93.606 


1.0 


7.624 


34.757 


119.714 



Table 2: Dependence of the instanton density n+ on the coupling parameter In the table 
the density n+ is given for three different values of e'^ {L = 1). One should note that the 
defect a has quite a large error, therefore the resulting values for the density at /5i = are 
different from the one in Tab. ^ due to this statistical error. A is the intercept from the 
combined fit of the free energy to a form F = A + N[ln{N/ L'^) — 1] + A^lnn + a In A^! and 
it is a small quantity, which can be neglected in the further discussion. 





a 


b 


0.1 

0.01 
0.001 


1.94273 
3.29677 
4.29687 


0.119835 
0.331814 
0.607095 



Table 3: Fit of the equilibrium density to the form given by the Debye-Hiickel approximation. 
The table shows the results for the linear fit to the form lnn+ = a + 61n/(/3i), where 
ln/(/3i) = |[(1 + (3i) ln(l + pi) + (1 - /3i) ln(l - /3i)]. The increase of the intercept when e'^ 
becomes smaller is due to the fact the the density for A^c = 2 is divergent. All values are for 
/3 = 1. 
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Figure 1: The 'dilute' regime for Nc = 3: Concentrating around the center-of-mass points Xi, 
X2, and 0:3 zindons of three different colors are closely grouped together to form instantons 
with large distances Rij between them. The configuration corresponds to a dilute gas of 
non-overlapping instantons. 
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Figure 2: Stability of the Metropolis algorithm for = 2 [e'^ = 0.001] (on the left) and 
A'^c = 3 (on the right, bold line) A^^ = 4 (on the right, thin line). The graphs show the free 
energy improved step by step by the Metropolis algorithm until a plateau of 'important' 
configurations is reached. In all three runs we used iV_|_ = 30 and L = 10. In all cases the 
interaction between instantons and anti-instantons is switched off (/3i = 0). 
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Figure 3: Snapshots of zindon configurations at the plateau after 500.000 Metropohs steps: 
The figure shows the spatial distribution of zindons for = 2 (on the top ), Nc = 3 
(bottom left) and iVc = 4 (bottom right). The filled and empty circles and squares symbolize 
different colors. Only for Nc — 2 a, condensation into neutral pairs is observed; for higher 
Nc the free energy is finite though many distinct 'color-neutral' A^c-plets are seen. Only 
instanton-zindons are shown. 
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Figure 4: Divergent/convergent behavior of the free energy: The free energy is plotted 
versus ln(ln(100/e)) for Nc = 2 (top), Nc = 3 (bottom left), and A^^^ = 4 (bottom right). 
One observes a divergent behavior for N,, = 2, while the free energy is finite for N^. = 3 and 
Nc — 4:. The last point in the plots with Nc — S and A'c = 4 is calculated using e = 0. 
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Figure 5: Density of the Coulomb gas {N^ = 2): The curves in the figure show the free 
energy F versus the number of instantons using 20.000 Metropohs steps until equilibrium 
is reached. From thin to bold we use the lengths L = 0.5, L = 0.6, L = 0.7, L = 0.8, 
and L = 0.9. The smooth solid curves show the combined fit with the scaling function 
F{L, N) = A^[ln(L^) — (In(A^) — 1 — ln(?7,+))] + a, with a single a and n+ given to fit all five 
curves simultaneously. If we fit each curve separately we get from thin to bold the following 
values for the density n+ = 62.341, = 64.526, = 66.027, n+ = 63.857, and n+ = 
63.836. In case of the combined fit we get for the density the value ?^+^bcst — 

64.117± 1.193, 

where the error is the standard deviation derived from the five separate measurements of n+ 
given above. The cutoff parameter in this figure is e'^ = 0.001. The defect a is 0.481. 
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Figure 6: Check of the Debye-Hiickel approximation to the Coulomb gas [N^ = 2). The 
sohd hues from thin to bold give the measured value for a Coulomb gas with /? = 0.1 and 
lengths L = 2, L = 3, L = 4, L = 5, L = 6. The thin smooth solid lines show the prediction 
of the scaling law, which is F(N,L) = (2 - p)N[ln{L'^) - (In(A^) - 1 - ln(n+))] + a, where 
in this case is the number of instantons only, as we do not consider any anti-instantons 
here. The four curves yield the value for n+ = 1.06787 ± 0.0024, which is in agreement with 
the prediction of the Debye-Hiickel approximation, which gives n+^oebye = 1.0370. 
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Figure 7: Fit of the equilibrium density to the form given by the Debye-Hiickel approxi- 
mation. According to the Debye-Hiickel approximation the logarithm of the density lnn_|_ 
is linear in ln/(/?i), where ln/(/?i) = + /3i) ln(l + /3i) + (1 - A) ln(l - In the 

figure the data for lnn+ are shown versus ln/(/?i), where the triangles belong to the sample 
with e'^ = 0.1, the squares to the sample with e'^ = 0.01 and the circles to the sample with 
e'^ = 0.001. The solid lines are linear fits. The corresponding fit constants are given in 
Tab. I 
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Figure 8: Free energy versus the number of instantons (equal to the number of anti- 
instantons), N/2 ~ N_ — N+. The figures show the free energy versus N/2 for N^ = 3 
(on the left) and Nc — 4: (on the right). The solid hues from thin to bold are for Pi — 0.0, 
Pi — 0.2, and Pi — OA. The dashed hne is the free energy for the value Pi — 1, which shows 
all signs of a collapse. 




Figure 9: Prediction of the Debye-Hiickel approximation for the density at /? = 1. The solid 
hne is the Debye-Hiickel prediction for Nc — 2, the dotted hne for A^^. = 3, the dashed-dotted 
hne for N^ — 4:, and the dashed line shows the limit Nc ^ oo. The lines end at critical 
couplings. 



35 




0.0 0.3 0.6 0.9 1.2 1.5 



p/<R> 




Figure 10: Instanton size distribution following from random (non-interacting) zindons are 
displayed for Nc — 2 (on the top), Nc — 3 (bottom left), and Nc — 4: (bottom right) 
for = N_ — 8. The histograms show the size distributions 'geometrically' defined by 
grouping zindons according to their dispersion p. The stars show the result for the 'lattice' 
method using a 100 x 100 grid. 
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Figure 11: Instanton size distributions with zindon interactions switched on are displayed 
for Nc = 3 (left), and A^^^ = 4 (right) for A^_|_ = N_ = 8. The histograms show the 'geometric' 
size distributions, and the solid hues the size distributions seen by the 'lattice' method, again 
using a 100 x 100 lattice. The instanton-anti-instanton coupling constant for the solid lines 
and the histograms is /?i = 0.0. The stars and the dashed lines show the 'geometric' and 
'lattice' size distributions for the case of f3i = 0.5. Instanton sizes are plotted in units of the 
average separation (R). The arrows show the maxima of the corresponding 'geometric' size 
distributions obtained in the non-interacting case of Fig. 
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Figure 12: Topology of a sample with 3 instantons and 3 anti-instantons {Nc = 3) seen by 
the 'lattice' method. The sample has been generated with the full partition function at the 
plateau for Pi — 0.5. On the left we show the the zindon content of the configuration. The 
triangles connect zindons of three different colors forming instantons (solid lines) and anti- 
instantons (dashed lines), as found by the 'geometric' method. Filled symbols correspond 
to zindons, open symbols to anti-zindons. The diamonds, squares and triangles denote the 
different colors. The stars show the centers and the circles show instanton sizes, as found by 
the 'lattice' method using a 100 x 100 grid. The method has been applied to the topological 
charge density displayed as a 3-dimensional plot on the right. 
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